Fermion mass and the pressure of dense matter 
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Abstract. We consider a simple toy model to study the effects of finite fermion masses on the 
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The role of finite quark masses in QCD thermodynamics has received increasing 
attention in the last few years. In the case of cold and dense QCD, it was generally 
believed that effects of nonzero quark masses on the equation of state were of the order 
of 5%, thereby yielding only minor corrections to the mass-radius diagram of compact 
stars yj]. In fact, mass, as well as color superconductivity gap, contributions to the 
pressure are supressed by two powers of the chemical potential as compared to zero- 
mass interacting quark gas terms. Therefore, assuming a critical chemical potential for 
the chiral transition of the order of a few hundred MeV, naively those terms should 
not matter. However, recent results for the thermodynamic potential to two loops have 
shown that corrections are sizable, and may dramatically affect the structure of compact 
stars [0]. Moreover, the situation in which mass (as well as gap) effects are significant 
corresponds to the critical region for chiral symmetry breakdown in the phase diagram 
of QCD. Hence, not only the value of the critical chemical potential will be affected, 
but also the nature of the chiral transition. In particular, if the latter is strongly first-order 
there might be a new class of compact stars, smaller and denser, with a deconfined quark 
matter core |3Q. Of course, contributions due to color superconductivity |4|] as well as 
chiral condensation [5] will also affect this picture. 

In what follows, we study a simple toy model - cold and dense Yukawa theory - to 
investigate the influence of fermion masses on the pressure. Here, we present a two-loop 
calculation of the pressure with massive fermions in the modified minimal subtraction 
(MS) renormalization scheme |6|], and briefly comment on possible implications to 
the physics of condensates in the core of neutron stars and effective models for color 
superconductivity. Higher-order corrections and a thorough analysis of renormalization 
group effects will be presented elsewhere [7|]. 

We consider a gas of massive fermions whose interaction is mediated by a real scalar 
field, 0, with an interaction Lagrangian of the Yukawa form, Jfj = g ty/ty/^ , where g is 
the coupling constant. In the zero-temperature limit, the perturbative pressure results in 



a power series of ay = g 2 /An. Q Up to 0((Xy), the first non-trivial contributions to the 
pressure are given by the free massive gas term, Pq, and the "exchange diagram", P\ . 
Using standard methods of field theory at finite temperature and density [8], one can 
derive the free gas pressure for fermions of mass m, obtaining in the zero-temperature 
limit the following form: 
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where jl is the chemical potential and Pf = V V 2 — m 2 denotes the Fermi momentum. 
The O(cty) renormalized correction reads [@]: 
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where w = — m 2 ln[(/j. + pf)/m] and A is the renormalization scale in the MS 
scheme. 

Fig. 1 illustrates the effect of modifying the mass on the total pressure to O(ay), 
P = Pq + Pi. The choice of range for /i, and accordingly for the masses, are inspired 
by the scales found in the case of QCD [2]. In the same vein, the coupling is fixed to 
OCy = 0.3. It is clear from the figure that mass corrections bring significant changes to the 
pressure, even in the absence of renormalization group (RG) running for the coupling 
and the mass. The figure also shows the dependence on the renormalization scale A. The 
values chosen are motivated by the ones which appear in QCD, as before. Although the 
effects of varying A appear to be relatively small, it would be premature to conclude that 
this feature will remain after implementing the RG flow. In fact, the results presented in 
Fig. 2 most probably underestimate the scale dependence of the full correction, since not 
only the coupling but also the mass will run with A. In the Yukawa theory, in contrast 
to QCD, the effect will become larger as we increase the chemical potential. For fixed 
coupling, larger values of A yield larger modifications in the pressure. However, after the 
inclusion of RG running, this behavior can be mantained, as should be the case here, or 
become the opposite, as is the case in QCD, depending on the sign of the beta function. 
Since the A-dependence comes from the term ~ m 2 ayln(A/m) in ©, there will be a 
competition between the behavior of the renormalization scale A and that of m and ay 
as functions of /i. 

Even at two loop order mass effects bring into play logarithmic corrections originated 
in the MS subtraction scheme. As usual, they bring about a non-physical dependence on 
the renormalization scale A, since one has to cut the perturbative series at some order. 
Higher-order computations in this framework are in progress D7|], and will give a better 
handle on the choice of this scale, which in our case should be a function of /i and m. 
On the other hand, one can also choose the scale in a phenomenological way in a given 
model, imposing physical constraints to the equation of state, as was done in Ref. PO to 
model the non-ideality of QCD at finite density with massless quarks. 



Since we are concerned only with the zero-temperature limit, there are no odd powers of g coming from 
resummed contributions of the zero Matsubara mode for bosons in the perturbative series. 
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FIGURE 1. Pressure normalized by the free fermion gas pressure as a function of the fermion chemical 
potential. Left: A = 2/1 and different values of the fermion mass. Right: m = 100 MeV and different values 
of the renormalization scale A. 



The points discussed above might be relevant in the study of effective models for 
the cold and dense matter found in the interior of compact stars, especially because 
the effects seem to be significant near the critical region. In the context of the NJL 
model, e.g., it was shown that a self-consistent treatment of quark masses strongly 
affects the competition between different phases [5]. And the mechanism of pairing 
in color superconductivity will certainly be influenced [@] by the running of nonzero 
quark masses. The investigation of these issues, as well as the effect of nonzero fermion 
masses in the formation of other condensates in neutron star matter, is under way [|7|]. 
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